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Abstract 

In this paper we are going to describe the results of the computer ex- 
periment, which in principle can rule out the Riemann Hypothesis. We 
use the sequence Ck appearing in the Baez-Duarte criterion for the RH. 
Namely we calculate Ciooooo with thousand digits of accuracy using two 
different formulas for with the aim to disproof the Riemann Hypoth- 
esis in the case these two numbers will differ. We found the discrep- 
ancy only on the 996 decimal place (accuracy of 10^^^^). The reported 
here experiment can be of interest for developers of Mathematica and 
PARI/GP. 



1 Introduction 

With the advent of the computers era the computing machines have been used to 
prove mathematical theorems. The most spectacular examples of such a use of 
computers were proofs of the four color theorem [2], [3] and of the Kepler conjecture 
about sphere packing in three-dimensional Euclidean space [12]. It seems to be 
not possible to use computers for the proof of the Riemann Hypothesis (RH), but 
refutation of it by numerical calculations seems to be plausible. 
The Riemann Hypothesis says that the series 

^(^) = E;^' is = ^ + ^t, ^{s}>i) (1) 

n=l 

analytically continued to the complex plane in addition to trivial zeros ({—2n) = 
has nontrivial zeros C{pi) = in the critical strip < 3ft{s} < 1 only on the critical 
line: ^{pi} = \ i.e. pi = ^ + i'ji, see e.g. the modern guide to the RH [29]. In the 
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same book there is a review of failed attempts to prove RH in Chapter 8. Presently 
the requirement that the nontrivial zeros are simple C'{pi) 7^ is often added. 

The first use of computers in connection with RH was checking by Allan Turing, 
whether the nontrivial zeros of ({s) have indeed real part ^ [3l]. Turing was con- 
vinced that RH is false. Let us quote the sentence from the first page of his paper: 
"The calculation were done in an optimistic hope that a zero would be found off the 
critical line" , but up to t = 1540 Turing has found that all zeros are on the critical 
line. The present record belongs to Xavier Gourdon [llj who has checked that all 
10^^ first zeros of the Riemann ({s) lie on the critical line. Andrew Odlyzko checked 
that RH is true in different intervals around 10^° [23], lO^^ [2l], lO^^ [27], but his 
aim was not verifying the RH but rather providing evidence for conjectures that 
relate nontrivial zeros of C(s) to eigenvalues of random matrices. In fact Odlyzko 
has expressed the view, that off critical line zeros could be encountered at least at t 
of the order 10^*^ , see ^, p. 358]. Asked by Derbyshire "What do you think about 
this darn Hypothesis? Is it true, or not?" Odlyzko replied: "Either it's true, or else 
it isn't". Also other famous mathematicians John E. Littlewood and Paul Turan 
have not believed RH is true. 

There were several attempts to use computers to disprove some conjectures re- 
lated to RH in the past. Sometimes it was sufficient to find a counterexample to 
the given hypothesis, sometimes the disprove was not direct. For example, the 
Haselgrove pTS] disproved the Polya's Conjecture stating that the function 

L(x):=5^A(n) (2) 

n<x 

satisfies L{x) < for s > 2, where X{n) is Liouville's function defined by 

A(n) = (-l)'-W 

where r(n) is the number of, not necessarily distinct, prime factors in n = ■ ■ ■ p^"^^-), 
with multiple factors counted multiply: r = ri -|- . . . + r„. From the truth of the 
Polya Conjectures the RH follows, but not the other way around. The Haselgrove 
proof was indirect, and in 1960 Lehman |16j has found on the computer explicit 
counterexample: L(906180359) = 1. 

The next example is provided by the Mertens conjecture. Let M{x) denote the 
Mertens function defined by 

M(x) = 5^/iH, (3) 

n<x 

where fi{n) is the Mobius function 

for n = 1 
jj,{n) = { when p^ln 




From 



when n = piP2 . . .pr 
|M(a;)|<x^ (4) 



Failed attempt to disproof the Riemann Hypothesis 



3 



again the RH would follow. However in 1985 A. Odlyzko and H. te Riele 
disproved the Mertens conjecture, again not directly, but later it was shown by J. 
Pintz [30] that the first counterexample appears below exp{3.21 x 10^^). The upper 
bound has since been lowered to exp{1.59 x 10"^°) [H]. 

Especially interesting is the value of the de Bruijn-Newman constant A, see e.g. 
§2.32 (pp. 203-205) in [10]. Unconditionally A < 1/2 and the Riemann Hypothesis 
is equivalent to the inequality A < 0. The fascinating run for the best lower bound 
on A ended with the value A > — 2.7x 10^^ obtained by Odlyzko [26]. Such a narrow 
gap for values of A being compatible with RH allowed Odlyzko to make the remark: 
"the Riemann Hypothesis, if true, is just barely true" . 

In 1997 Xian-Jin Li proved [T7j|, that Riemann Hypothesis is true iff the sequence: 

where 

as) = ls{s-l)T (0C(3) 

fulfills: 

A„ > for n = 1,2,... (5) 
The explicit expression has the form: 

An = $^(l-(1-1/Pr). (6) 

p 

K. Maslanka [21], [20] has performed extensive computer calculations of these con- 
stants confirming ([s]). 

Let us mention also the elementary Lagarias criterion [15]: to disprove the RH 
it suffices to find one n that has so many divisors, that: 

J2d>H,, + exp{H^) \og{H^). (7) 

d\n 

The Lagarias criterion is not well suited for computer verification (it is not an easy 
task to calculate Hn for n ~ 10^"°°'^° with sufficient accuracy) and in [7] Keith Briggs 
has undertaken instead the verification of the Robin [31j criterion for RH: 

RH < e'^nloglog(n) for n > 5040 (8) 

d\n 

For appropriately chosen n Briggs obtained for the difference between r.h.s. and 
l.h.s. of the above inequality value as small as e~^^ ~ 2.2 x 10"^, hence again RH is 
very close to be violated. 

In this paper we are going to propose a method which in principle can provide a 
refutation of the RH. The idea is to calculate some number with very high accuracy 
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(one thousand digits) in two ways: one without any knowledge on the zeros of ({s) 
and second using the exphcit formula involving all pi. Despite some estimation 
presented in Sect. 3 indicating that the discrepancy could be found merely at much 
higher than a thousand decimal place we performed the calculations in an optimistic 
hope that we will find the discrepancy between these two numbers, paraphrasing 
the sentence of Turing. There is a lot of number theoretic functions defined often in 
an elementary way being expressed also by the "explicit" formulas in terms of zeros 
of the ({s) function. Let us mention here the Chebyshev function 

ij{x) = Mn) (= log(lcm(2, 3, ■ ■ ■ [x\ ))), 

n<x 

where von Mangoldt function A{n) is defined as 

A{n) = 

The explicit formula reads, see eg. [33J: 



logp for n = p" 
in other cases 



C'(0) ^ xP^ , / ^ 1 , 1 \ a;/", , 

all zeros pi ' nontr. zeros pi 

Also the Mertens's function has the explicit representation (last term is compris- 
ing contribution from all trivial zeros) |35j : 



n<x nontr. zeros pi \i / n=l ^ ' 



2n ^ 

{2ny.nC{2n + 1) 



The problem with these series is they are extremely slow to converge because the 
partial sums oscillate with amplitudes diminishing at very slow rates. For example 
?/;(1000001) = 999586.597..., while from iQ summing over 5,549,728 zeros gives 
999587.15 . . ., thus relative error is 0.000055. 

In the computer experiment reported here we were able to get discrepancy less 
than 10~^^^ between quantity calculated from generic formula and from explicit one 
summed over only 2600 nontrivial zeros computed with 1000 significant digits. 

2 The Baez-Duarte criterion for the Riemann Hypothesis 

We begin by recalling the following representation of the (^{s) function valid on the 
whole complex plane without s = 1 found by Krzysztof Maslanka [19] : 

.-1^ r(i-f) « 
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where 



j=0 ^"'^ j=0 ^"'^ 



TT 



2j+2 



5 



2i+2- 



(12) 



The expansion (11) provides an example of the analytical continuation of (jl]) to 
the whole complex plane except s = 1. Since Ak tend to zero fast as A; — > oo the 
expansion (11) converges uniformly on the whole complex plane Based on the 
representation (11) Luis Baez-Duarte in has proved that RH is equivalent to the 
statement that 

Ck = , V£ > (13) 



where 



If additionally 



Ck 



j=0 



j7C(2j + 2)' 



(14) 



Ck = {k-'") (15) 

then all zeros of C(s) are simple. Baez-Duarte has shown unconditionally (regardless 
of validity of the RH) slower decrease Ck = 0{k~^/'^). The plot of Ck for k = 
1,2,..., 100000 is presented in the Fig.l. 

The explicit formula expressing the Baez-Duarte sequence Ck directly in terms 
of the zeros of ({s) can be written as a sum of two parts: quickly decreasing with 
k trend Ck arising from trivial zeros of ({s) and oscillations Ck involving complex 
nontrivial zeros: 

Cfc = Cfc + Ck (16) 
Because the derivatives ('{—2n) at trivial zeros are known analytically: 



C'(-2n) 



-l)"C(2n + l)(2n)! 



(17) 



Maslanka in [IB] was able to give the closed expression for trend stemming from 
zeros pn = —2n: 



Ck 



(2vr) 



r(fc + i)r(m) 



-l)'"(27r) 



2m 



T{k + m + l)T{2m-l) C(2m - 1) 



Baez-Duarte is skipping the trend Ck remarking only that it is of the order o(l) 
(Remark 1.6 in [5J). It is an easy calculation to show (see |36j), that for large k 



Ck 



(2vr)!l 
2C(3) P 



+ 0{k' 



(19) 



thus the dependence Ck = 0{k ^/^) in (13) is linked to the oscillating part c^. 
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-4.0x10 " - 



-8.0x10 




200000 



400000 



600000 



800000 



1000000 



Fig.l The plot of the Baez-Duarte sequence Ck for k G (1, 10^). The equation for 



envelope was obtained from the explicit equation (27): for large k the oscillating 



part Ck is dominant and for k > 100000 Cfc fits well between the red lines, for 

details see 



For Ck Baez-Duarte gives the formula [3]: 



pi 

where 



^^(^) - n (l - 3 (21) 

r=l 

is the Pochhammer symbol. Assuming zeros of C(s) are simple we can write: 

~''~' = 2k^C'{p,)Pk{p^/2y (22) 



An appropriate order of summation over nontrivial zeros is assumed in (20) and 
(22), see [Sj Theorem 1.5]. Because 

PAs) - (-^)'^(^) (23) 



Failed attempt to disproof the Riemann Hypothesis 



7 



collecting in pairs pi and p we can convert (22) to the form: 

k-1] 



Ck 



^ r(fc + i)r(f 



(24) 



^pi, S(p,)>0 

We have found that the numerical calculation of Pk{pi/2) in PARl/GP directly from 
above product (21) is much slower than the use of the T{z) functions (23). 
Baez-Duarte proves in jH Lemma 2.2] that 



lim Pk{s)k' 

k^oo 



1 



rn-s) 



(25) 



thus for large k we can replace + 1 by A; and transform (|22|) in the following way: 

kP'/^ 



Ck 



2k 



E77 



pi 



C{pi)Pk{pi/2)kP'/' 




CiPi) 



(26) 



Now we assume the RH: pi = Then we get for Ck the overall factor k 



the dependence following from RH, see (|13|) — multiplied by oscillating terms: 

1 



Ck 



^3/4 



\-yt>o 



(27) 



Using the formula (6.1.45) from yj: 

1 



lim 



\y\-*oo v27r 

assuming RH we obtain for large 7^ > 0: 



r(x + i?/)|e^^l^l|?/ 



r 



4^2^' 



27re 



-7r7i/4 fJi\ 4 



(28) 



(29) 



hence we get exponential decrease of summands in the sum (27) over nontrivial zeros 
giving Ck and (27) is very fast convergent. Because of that if RH is true the sum (27) 
will be dominated by first zero 71 = 14.13472514 . . . what leads to the approximate 
expression (for details see [36]): 



Ck 



A 



sin - -7i log(A;) A = 7.775062 ... x 10' 



2.592433 . 



(30) 



For large k the above formula (30) gives very fast method of calculating quite accu- 
rate values of Ck, orders faster than (14). 

In the following we will denote by the values of the Baez-Duarte sequence 
obtained from the generic formula ( 14 ) and by the values obtained from explicit 
formula (16), i.e. in fact from (18) and (27) as no one zero off critical line is known. 



8 



Marek Wolf 



3 The scenario of violation of the Riemann Hypothesis 



The condition (15) means that the combination k^Ck should be contained between 



two paralel hnes ±C for all k, where C is the constant hidden in hig-O in (15). 



The violation of the RH would manifest as an increase of the amplitude of the 

3 3 

combination k^Ck and for sufficiently large k (depending on C) the product k^Ck will 
escape outside the strip ±C (if RH is true we can take C = A = 7.7751 . . . x 10^^). 



We will discuss below the case (15) of simple zeros {('{s) ^ 0) requiring Ck = 
0(fc-3/4). 

The derivatives C{pi) the denominator of (26) does not pose any threat to 
RH. First of all it does not depend on k, thus hypothetical extremely small val- 
ues of C{pi) will oiily change the constant hidden in big-C in (15). Second, this 



derivative is taking moderate values for zeros used by us: the smallest ^ 
0.032050162. . . at 71310 
71773 = 2275.06866478. 



2 +7;) was 

= 1771.212945 ... and the largest was 7.7852581838 ... at 
We remark that Lehman in [TB] makes the conjecture: 



1 



C'iPi) 



0{p. 



where < z/ < 1. 



(31) 



Some rigorous theorems about the possible large and small values of C'(Pz) proved 
under the assumption of RH can be found in We have checked, that for first 
5549728 nontrivial zeros of C'(pO largest derivative was 9.38127677 at 75376610 
and the smallest was 0.001028760514 at 74161179- 

Let us suppose there are some zeros of C,{s) off critical line. Next let us assume 



that we can split the sum over zeros pi in (26) in two parts: one over zeros on critical 
line and second over zeros off critical line. This second sum should violate the overall 
term present in the first sum. Let p[°'' denote the zeros lying off critical line 

("o" stands for "off"): pf^ = \ -^5i + i'')\°\ 5i > Q (as it is well known the nontrivial 
zeros are symmetric with respect to the critical line zero hence the combination ±5; 

plus there are appropriate conjugate zeros below real axis). In the factor |r(l 

1 

the off-critical line zeros will lead only to the change of 7^* 



7;^^'^'^^ in (|29D and in 



k- 



PL" 



view of exponential decrease present in (29) the possible violation of RH in (26) will 
manifest through the terms k^'^"^'^ 



The combination 1/ki 



+S,/2 



poses no problem as it leads to faster than required in (|15|) decrease of some of 

3 



terms in the sum for c^. But the combination l/fcj"*^'/^ leads to violation of ([l5| 
and we want to elucidate how it can be detected. Say we want to compare with 
df, with accuracy e, where we are interested in values of e of the order e = 10~^°'"°. 
The expression for is a finite sum and we can in principle calculate its value in 
PARI with practically arbitrary exactness (however for really large k it can take 
years of CPU time). Although contains information coming from all zeros, to see 
infiuence of the first off critical line zero the value of sufficiently large k has to be 
examined. The sum for is infinite and we expect that to get accuracy e we have 



to sum in (27) up to Z = L given by (as we skip |7i| ' we will skip also 7r/4 as our 
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consideration are not rigorous in general): 

e^^e""^^, hence 7^ f« — log(e). (32) 

Because values of the imaginary parts 7^^°'' of the hypothetical zeros off critical line 
should be extremely large, perhaps even as large as 10^° , we suppose that ■jl < 
'-)\°\ The contribution of 7^°'' is present in c^, but will not be present in the explicit 
sum for cut at L. To detect discrepancy between and larger than assumed 
accuracy e sufficiently large value of = i^' is needed. The point is that fc^^^c^ will 
escape outside the strip ±C for sufficiently large k = K and the value of such K we 
can estimate analyzing the explicit formula for c^. 

We can estimate value of the index K from the requirement that the term K^'-^'^ 
(i^"^/^ is present in front of the sum for Ck) will defeat the smallness of the term 
r(3/4 — '-fl"^ /2) and together their product will overcame the first summand in 
(26) corresponding to 71. In other words in the series (26) all terms up to 'yi°^ 



monotonically and fast decrease but the terms corresponding to zeros off critical 
line can be made arbitrary large for sufficiently large k. The condition for such a K 
is roughly: 

K^i/'^e-^^^"^ > C thus > CV'^'"'/^'. (33) 

Because 61 can be arbitrarily close to zero and, as we expect, 7!°'' is very large the 
value of K will be extremely huge — larger than famous Skewes number and will look 
something like 10^" . Thus it is practically impossible to disproof RH by comparing 



and c^. Maslanka has given in [TH] discussion of possible violation of (15) and he 



also came to the pessimistic conclusion that disproving RH by comparing and is 
"far beyond any numerical capabilities", see pp. 7-8. We wanted to find agreement 
between cfooooo c^ooooo within one thousand digits and to our surprise the first 
attempt to calculate cfooooo resulted in the difference already on the 87 place. We 
started to struggle with numerical problems to improve the accuracy and finally we 
got 996 digits of cfooooo '^100000 same. 



4 The computer experiment 

The idea of the experiment is to calculate to high precision the values of cfooooo 
cfooooo fi^*^ discrepancy between them. We calculated one value cfooooo 



from the generic formula (14), which contains contribution from all zeros of C{s), 
even hypothetical zeros with 3?(pz) 7^ |. Because C(2'^) very quickly tend 1 to get 
the firm value of Ck it is necessary to perform calculations with many digits accuracy. 
Additional problem is fast growing of binomial symbols. We performed calculation 
of cfooooo using the free package PARI/GP |32]- This package allows to perform very 
fast computations practically of arbitrary precision. We have set precision to 100000 



decimals and below in Table I are the partial sums of ( 14 ) recorded after summation 



of 10000, 20000, ... 100001 terms. In the middle of computations the partial sums 
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for c?ooooo were of the order lO^oioo to drop finally to 1.609757993 . . . x 10"^ after 
summing up all 100001 terms. Separately we have repeated this calculations with 
precision set to 150000 places and we found the difference only from the 69900 place. 



Because the sum (14) is finite we got accurate say 50000 digits from this generic 
formula for cfQQQQg. This was an easy part, which took approximately 14 hours for 
precision 100000 decimals and almost 20 hours for precision 150000 digits on the 
AMD Opteron 2.6 GHz 64 bits processor. 



Table I 



n 


v^n / -,N7/1UUUUU\ 1 
2^j=0y ^) \ J Jc(2j+2) 


10000 


5.65168726144550 x lO^^^^^ 


20000 


4.00927204946289 x lO^^^^g 


30000 


6.08771775660005 x 10^6526 


40000 


5.17938759373151 x 10^9225 


50000 


1.26030418446100 x 10^^^^^ 


60000 


3.45292506248767 x 102^225 


70000 


2.60902189568574 x 10^6526 


80000 


1.00231801236572 x 10^1^29 


90000 


6.27965251271723 x 10^^^^^ 


100000 


1.60975799392038 x 10"^ 



Now we turn to the calculation of c^qqqoo from the explicit formulas (18) and (27) 



which are infinite sums. It is important not to make the replacement fc + 1 — >■ /c even 
for k = 100000 if we want to get accuracy of the order of 1000 digits. The series in 



(18) decreases very fast to zero and it is very easy to get arbitrary number of digits 
of Ck using the PARI/GP procedure sumalt. 

Next we want to calculate Ck from (27) with 1000 digits accuracy. From e~'^""'^^ = 
]^Q-iooo 2931.7 and a glance at the list of zeros of ({s) (e.g. [25]) gives 

/ = 2402 because 72402 = 2931.0691 . . .. First 100 zeros of ({s) accurate to over 1000 
decimal places we have taken from Next 2500 zeros of ({s) and derivatives C'{pi) 
with precision 1000 digits we decided to calculate using the built in Mathematica 
V.7 procedures ZetaZero [m] and numerical differentiation ND[. . .]. After a few 
days we got the values of 7; and Ci.Pi)- We have checked using PARI/GP that these 
zeros were accurate within at least 996 places in the sense that in the worst case 
|C(p«)l < 10"^^^ / = 1,2, . . . ,2600. The formulas ^ and (fl8| were implemented 



in PARI with precision set to 1000 digits and to our surprise we obtained that the 



values of c\qqqqq and c^qqqoo coincide only up to 87 place: 



10-^6 > 



-100000 



-100000 



> 10 



-87 



(34) 
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The arguments given earlier in Sect. 3 suggested that no discrepancy should be 
found between and c| for k = 100000 hence we have paid attention to the numeri- 
cal inaccuracies in the computation of the derivative C{pi) ^ possible explanation. 
We have played with different options 20, 30, 40 terms and WorkingPrecision in 
ND [ . . . ] , but finally we got only the moderate improvement: the difference between 



•^100000 '^100000 shifted to 105 place: 



10 



'105 



< 



-^100000 



-^100000 



< 10 



'104 



(35) 



Because in (27) gamma functions can be calculated in PARI with practically 



arbitrary digits of accuracy the only way to improve accuracy of calculation of 
^100000 is to find a reliable method of calculating C'(s) with certainty that say all 
1000 digits are correct. From ([T]) it is easy to obtain the modified expression for 
zeta: 

n=l 



1 



This series is absolutely convergent for 3f?(s) 
by term: 



a > and can be differentiated term 



C'(^) 



log(2)2i- 



n-s\2 



E 

n=l 



-1 



,n-l 



+ 



1 



n=l 



[-l)"log(n) 



(37) 



PARI contains the numerical routine sumalt for summing infinite alternating 
series in which extremely efficient algorithm of Cohen, Villegas and Zagier [8] is 
implemented. As this authors points out on p. 6 their algorithm works even for 
series like (36) with s complex — (36) is alternating only when s G M. We have 



used this routine sumalt outside scope of its applicability with success to calculate 



C'ipi) from (37) with precision set to 1000 digits and then calculated c^ooooo from 
^ and (ITS) 



The result was astonishing: the difference between c^ooooo ^^"^ ^looooo 



appeared on the 625 place: 



-100000 



1 



-100000 



< 10 



-624 



(3^ 



Because we expect possible violation of RH should manifest at much larger k we 
were looking for the way to still improve the accuracy of C'iPi)- We decided to make 
a frenzy think: we calculated again C'{pi) using sumalt with zeros having 1000 digits 
but this time with precision set to 2000 (however values of C{Pi) were stored only 
with 1000 digits). Thus the aim was to enlarge the number of terms summed in 



(37), or in fact the number of iteration performed inside sumalt until the prescribed 



accuracy is attained. After 18 hours on AMD Opteron 2.6 GHz we got the results. 



And now bingo! The first 996 digits of cfgoooo c^ooooo where the same: 



10 



-996 



< 



-100000 



-100000 



< 10 



-995 



(39) 
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0.01 



y(n) - 



0.001 — 








200 400 600 800 1000 



1 \ \ \ r 

1200 1400 1600 1800 2000 2200 2400 2600 

n 

Fig. 2 The plot of y{n) for c^qoooo obtained by three methods. The plot (a) is the 
best result obtained with Mathematica with the procedure ND [ . . . ] with the option 
terms=30 and WorkingPrecision=1000. The accuracy was 105 digits and the 
curve (a) departs from the red line at the zero 7149 = 317.73480594237. . . for which 
exp(-7r7i49/4) = 4.193107483 x 10"^°^ The plot (6) was obtained with PARI from 
the formula (37) using the procedure sumalt with precision set to 1000 digits. 
The y{n) reaches plateau 0.0006856926750 ... at zero 71412 = 1884.00577834967 . . ., 
for which exp(— 7171412/4) = 2.391868726 x 10^^^^. The curve (c) is the same 
as (6) but derivatives + ^li) were calculated with sumalt and precision set 
to 2000 digits. Curve (c) reaches the plateau 0.00043191361 ... at the zero 
72430 = 2960.033617812 ... for which exp(- 7172430/4) = 2.298783954 x 10-^°^°. 



Failed attempt to disproof the Riemann Hypothesis 



13 



Because we got the precision (39) it is a posteriori proof that 1000 digits of 



derivatives C'{pi) were calculated correctly from (37) by the PARI procedure sumalt 
with precision set to 2000 digits. 

In the Fig. 2 we present summary of these computer calculations. Since it is not 
possible to plot using standard plotting software as small values as 10~^°° on the 
y-axis we present in the Fig. 2 the following quantity measuring the distance from 



cf 00000 ^he partial sums over zeros in (27) and decreasing with number of zeros 



included in the sum: 



y{n) = log 



^3/4 



, 1=1 



- ct 



(40) 



where k = 100000 and the absolute value is necessary as the differences between suc- 
cessive approximants to c^ooooo "^looooo changes sign erratically. The consecutive 
terms in the series (27) behave like e~'^""-^^ hence we expect that y{n), by analogy 



with well known property of alternating series with decreasing terms, should behave 
like the first discarded term: 

4 

y{n) . (41) 

The Fig.2 confirms these considerations: initially y{n) for all curves follows the 



prediction (41) and starting with n for which the values of the C'(Pn) ^ire incorrect 



adding further terms does not improve accuracy. Because of the exponential decrease 
of r(| — y) the contribution of further terms is suppressed and horizontal lines in 
Fig.2 are determined by the first corresponding to the bad value of the derivative 

C'(Pn). 



5 Final remarks 



Although we have reached the agreement between c^qoooo '^looooo places: 
.9 _ 1 Rno7^7oo ^ 980 digits ... 2913696 30140 x 10"^ 



Ciooooo = 1-60975799, 



'100000 



1.60975799 . . . 



980 digits 



. . .291369629833 x 10" 



it means nothing about the validity of the RH. The refutation of the RH by computer 
methods seems to be as difficult as the analytical proof of its validity. As a possible 
precaution let us mention the paper "Strange Series and High Precision Fraud" 
written by Borwein brothers [B]. In this paper we found a few striking examples 
of the approximate equalities correct to many, many digits, which finally are not 
identities. The most fraudulent is the following 



E 

n=l 



ne 



^163/9 



= 1280632 



(42) 
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which is vahd up to accuracy at least io^^°'^'°'^0'°oo but is not an identity. 

At the end let us remark that the most accurate experiment in physics is the 
measurement of the ratio of the electric charge of the electron e~ to the charge of 
the proton which is known to be something like e~/e^ = — 1 ± 10^^", see [35] . 
Physicist believe that e~ = — exactly. 

Acknowledgment: I would like to thank Dr. Keith Briggs and Prof. Krzysztof 
Maslanka for comments and remarks. 
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